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Abstract: We have in earlier papers presented an extension of Mueller's dipole cascade 
l/"") | model, which includes sub-leading effects from energy conservation and running coupling 

as well as colour suppressed saturation effects from pomeron loops via a "dipole swing". 

The model was applied to describe the total and diffractive cross sections in pp and j*p 
f"*- ■ collisions, and also the elastic cross section in pp scattering. 

In this paper we extend the model to describe the corresponding quasi-elastic cross sections 

in 7*p, namely the exclusive production of vector mesons and deeply virtual compton 

scattering. Also for these reactions we find a good agrement with measured cross sections. 

In addition we obtain a reasonable description of the ^-dependence of the elastic pp and 



quasi-elastic 7*p cross sections. 
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1. Introduction 

We have in a series of papers [1-3] presented an extention of Mueller's dipole cascade 
model [4-6] implemented in a Monte Carlo program, which includes sub-leading effects 
from energy conservation and running coupling, as well as colour suppressed effects from 
pomeron loops via a dipole swing mechanism. It also includes a consistent treatment of non- 
perturbative confinement effects, which suppress dipoles with large transverse extension. 

The advantage of a cascade model formulated in transverse coordinate space is the 
possibility to include effects of multiple collisions and saturation in a straight forward way. 
While analytic results have mainly been presented for the asymptotic behavior of total and 
diffractive cross sections, Monte Carlo simulations facilitate studies of non-leading effects 
and more quantitative results. A simulation of Mueller's initial model was presented by 
Salam in ref. [7]. Although giving finite results for the total cross section, this leading log 
evolution suffers from divergences for small dipoles, which caused numerical problems with 
very large gluon multiplicities and prevented simulations at higher energies. One important 
result from this analysis was the very large fluctuations in the evolution [8]. As the ratio 
between the elastic and the total cross sections is determined by the fluctuations in the 
scattering process, this implies that less fluctuations is needed in the impact parameter 
dependence, to reproduce the experimental data. As a result we found in ref. [3] that 
including the fluctuations in the evolution implies that the impact parameter profile is not 
as "black and white" as in analyses where only fluctuations in the impact parameter are 
taken into account. 

In the model described in refs. [1-3] we include a number of sub-leading effects, with 
the aim that we in the end will be able to describe not only the total and diffractive 
cross sections, but also to generate fully exclusive final states. The main ingredient in our 
model is energy conservation, which is included by assigning a transverse momentum to 
each emitted gluon given by the maximum inverse size of the neighboring dipoles. As a 
result this also implies that the singularities for small dipoles are avoided. Other features 
are saturation effects in the evolution through a dipole swing mechanism, and a consistent 
treatment of confinement and running coupling effects in both dipole emissions and dipole- 
dipole interactions. 

Taken together with a very simple model for the initial proton wavefunction, these 
features allow us to obtain a Lorentz- frame independent description of total cross sections, 
both for pp and DIS, using basically only two free parameters, a confinement scale r max and 
Aqcd [3]. The model gives a good description of measurements of the total and diffractive 
cross sections in pp and j*p collisions, and also for the elastic cross section in pp. In this 
paper we will continue our investigations with an analysis of exclusive production of vector 
mesons and real photons in 7*p. The aim is to further test our model, and in particular to 
study the effect of the fluctuations in the cascade. We also extend the analyses to include 
the ^-dependence of the (quasi-)elastic cross sections, including also elastic pp scattering, 
which in particular gives information about the properties of the incoming proton state. 

In the eikonal approximation the quasi-elastic 7*p collisions contain three elements: 
the virtual photon-dipole vertex, the dipole-proton scattering amplitude, and the vertex 



for the transition between the dipole and the final vector meson or real photon. Here the 
first component can be calculated perturbatively, although a hadronic component must be 
included at lower Q 2 -values. In an extensive study Forshaw et al. [9, 10] have analyzed the 
results obtained from a set of models for the dipole-proton scattering and for the vector 
meson wavefunctions, and compared them with experimental data. In this paper we want 
to carry out a similar analysis, but now use our dipole cascade model for the dipole-proton 
scattering. We are here particularly interested in effects of fluctuations in the cascade 
evolution, which are not included in the analyses by Forshaw et al. We also want to use 
this study to put constraints on the state of the incoming proton. 

We begin in section || with discussing the eikonal formalism for exclusive vector me- 
son production, whereafter we describe our model for dipole evolution and dipole-dipole 
scattering in section || and the models we use for the proton, photon and vector meson 
wavefunctions in section ||. In section || we retune the parameters of our model to data 
on total and elastic pp cross sections and the total 7*p cross section before we present our 
results on quasi-elastic rfp cross sections on section |6[ Finally we present our conclusions 
in section [7| 

2. Formalism 

2.1 The dipole cascade model and the eikonal approximation 

As discussed in the introduction, our model for pp collisions and DIS is an extension of 
Mueller's dipole cascade model [4-6]. In this formalism the probability per unit rapidity 
Y that a dipole (x, y) emits a gluon at transverse position z is given by 

dV a 2 {x-yf ... _ 3a s 

— — = — a z- -77- -77, with a = . (2.1) 

dY 2vr (x-zf(z-yf vr v ; 

The evolution of this cascade agrees with the leading order BFKL evolution. As a conse- 
quence, the total number of dipoles grows exponentially. This also implies a strong growth 
for the total cross section which, however, is tamed by taking multiple dipole interac- 
tions into account. The scattering probability between two elementary colour dipoles with 
coordinates (xj,y^) and (xj,yj) respectively, is given by 

a 2 [ / (x i -y j ) 2 (y i - Xj) 2 \] 2 
fij = f( Xi , yi \ Xj , yj ) = - pog ^- _ xi)2(yi _ ~y 2 j] . (2.2) 

Summing over all i and j this can give a large number, but the unitarized interaction 
probability determined by 

T(b) = l-e" E/ ^ = l-e" F (2.3) 

will never exceed 1. 

We note that the splitting probability in eq. ( |2.1| ) is singular for small dipole sizes 
x — z or z — y, but these small dipoles have a small probability to interact with the 
target, and the eikonal F = Y2 fij is finite. 



In the model developed in refs. [1-3], we extended Mueller's cascade model to include 
sub-leading effects from energy conservation and a running coupling, saturation effects not 
only in the dipole-dipole subcollisions but also within the individual cascades, and effects 
of confinement. These features are further discussed in sec. ||| 

The model is supplemented by a non-perturbative model for an initial proton in terms 
of three dipoles. In the eikonal approximation the total and the diffractive (including the 
elastic) cross sections are then given by 

a tot = 2Jd 2 b((l-e- F )}, a diS = Jd 2 b((l-e- F ) 2 }. (2.4) 

The diffractive cross section can be separated in elastic scattering and diffractive excitation: 

a el = fd 2 b ((1 - e~ F )) 2 , a diScxc = J d 2 b {((1 - e~ F ) 2 ) - ((1 - e- F )) 2 } . (2.5) 

Thus the separation between elastic and inelastic diffraction is determined by the fluctu- 
ations in the scattering amplitude. The average in eqs. ( |2.4D and (|2.5| ) is taken over the 
different incoming dipole configurations and different cascade evolutions, which thus give 
two separate sources for fluctuations. The differential cross section da e \/dt is obtained from 
the Fourier transform of the scattering amplitude 

2 

, with t = -q 2 . (2.6) 



fd 2 be icih {l-e- F ) 



dax __ J_ 
dt 4?r 

2.2 DVCS and exclusive vector meson production in ~f*p collisions 

We want to study the exclusive processes 

J*p^-Vp, V =j,p,ip,. .. (2.7) 

In the dipole model the virtual photon is split into a qq pair long before the collision. This 
dipole scatters elastically against the proton, and after the scattering the pair joins again 
forming a real photon or a vector meson. The formulation in the transverse coordinate plane 
makes it easier to study these pseudo-elastic reactions, and in the eikonal approximation 
the scattering amplitude is expressed in terms of three components: 

ZA x (s,b)=sYl I i'dzd 2 r*f h l(r,z)^ ih (r,z,Q 2 )a dp (s,r,b,z). (2.8) 

f,h,h 
Here r is the transverse size of the dipole, z and 1 — z the fractions of the photon or vector 
meson carried by the quark and antiquark respectively, and h and h their helicities. A 
denotes the photon or vector meson helicity, a dp is the dipole-proton scattering probability 
with b the impact parameter, and s the total energy squared. 

Neglecting the small contribution from the real part of the amplitude, the total cross 
section is given by 

<rx{l*P^Vp) = -Ljd 2 b\A x (s,b)\ 2 - (2-9) 

In analogy with eq. (^^) the differential cross section is obtained from the Fourier trans- 
form: 

-^=16^2 fd 2 be^ h A x (s,h) , witht = -q 2 . (2.10) 



3. The improved dipole cascade 

As discussed in the introduction, the model developed in refs. [1—3]. is an extension of 
Mueller's dipole cascade model, which includes sub-leading effects from energy conservation 
and a running coupling, saturation effects not only in the dipole-dipole subcollisions but 
also within the individual cascades, and effects of confinement. As mentioned above, an 
essential point is here that we include the effect of fluctuations in the dipole cascades in 
the calculation of the elastic or quasi-elastic cross sections. 

3.1 Non-leading perturbative effects 

3.1.1 Energy-momentum conservation 

It is known that the large NLO corrections to the BFKL evolution are reduced signifi- 
cantly if proper energy conservation is included in the leading order. In our model a small 
transverse extension is interpreted as a large transverse momentum. This interpretation 
is supported by the resulting analogies between the dipole chains in coordinate space and 
the chains in the LDC model, which is formulated in momentum space and interpolates 
smoothly between DGLAP and BFKL evolution. Taking energy-momentum conservation 
into account is most easily done in a Monte Carlo (MC) simulation. Conserving both 
light-cone components, p+ and p_, implies that we also satisfy the so called consistency 
constraint [11]. As small dipoles in our formalism correspond to large transverse momenta, 
energy conservation also gives a dynamical cutoff for the otherwise diverging number of 
small dipoles, and thus makes the MC simulation much more efficient. 

3.1.2 Running coupling 

In our simulations we also include non- leading effects from the running of a s , both in the 
dipole splitting and in the dipole-dipole scattering probability. In the dipole emissions the 
scale in the coupling is given by min(r, n, r2), where r is the size of the mother dipole which 
splits into n and ri- This is the most natural choice and is also consistent with recent 
NLO calculations [12-14]. For the dipole-dipole scattering the situation is somewhat more 
complicated with basically six different dipole sizes involved. We have chosen to use the 
scale min( \xi — y^ | , \xj — %fj \ , \xi — y ,• | , \y^ — Xj |) . In order to avoid divergencies the coupling 
is in all cases frozen so that a s (r) — > a s (min(r, r max )), where r max is the confining scale 



discussed in section 3.3 below. 



3.2 Saturation within the cascades 

Mueller's cascade includes saturation effects from multiple collisions in the Lorentz frame 
chosen for the calculation, but not saturation effects from gluon interaction within the 
individual cascades. The result is therefore dependent on the chosen Lorentz frame. In 
ref. [2] we improved our model by allowing (colour suppressed) recouplings of the dipole 
chain during the evolution, a "dipole swing" . The swing is a process in which two dipoles 
(aj^yj and (xj,yj) are replaced by two new dipoles (xi,yA and {xj,yj). The process 
can be interpreted in two ways. There is a probability 1/-/V,? that the two dipoles may 
have the same colour, and the quark at Xi and the antiquark at y^ form a colour singlet. 



In this case the best approximation of the quadrupole field ought to be obtained by the 
closest charge-anticharge combinations. Here the swing is therefore naturally suppressed 
by 1/N 2 , and it should be more likely to replace two given dipoles with two smaller ones. 
Secondly, we may see it as the result of a gluon exchange between the dipoles, which results 
in a change in the colour flow. In this case the swing would be proportional to a 2 , which 
again is formally suppressed by N 2 , compared to the splitting process in eq. ( |2. 1| ) , which 
is proportional to a = N c a s /ir. 

In the MC implementation each dipole is randomly given one of N 2 possible colour 
indices. Only dipoles with the same colour can swing, and the weight for a swing 
(xi,yi), (332,2/2) ~" * (^l j 2/2)5 ( a; 2,2/i) is determined by a factor proportional to 

(x 1 -y 1 ) 2 (x 2 -y 2 ) 2 , 31 . 



(a3i -y 2 ) 2 ( x 2-yi) 2 ' 

This implies that the swing favors the formation of smaller dipoles. The number of dipoles 
is not reduced by the swing, but the fact that smaller dipoles have smaller cross sections 
gives the desired suppression of the total cross section. Although not explicitely frame 
independent the results from the MC simulations are very nearly independent of the Lorentz 
frame used for the calculations. 

3.3 Confinement effects 

Mueller's dipole model is a purely perturbative process. It should therefore be applied to 
small dipoles, e.g. to heavy quarkonium states. When applying the dipole formalism to 
collisions with protons it is necessary to take confinement into account, in order to prevent 
the formation of very large dipoles. Confinement effects must also suppress long range 
interactions between colliding dipoles. In ref. [3] a consistent treatment of confinement 
was presented by replacing the Coulomb potentials in eqs. fl2.2| ) and ( |3.3| ) by screened 
potentials, with a screening length r max . 



Obviously the dipoles produced in the splitting process in eq. (2.1) cannot become too 
large, and it is natural to introduce a scale r max , so that larger dipoles are suppressed. In 
a similar way confinement must suppress long range interactions between colliding dipoles. 



The formula for fy in eq. (2.2) is just the two dimensional Coulomb potential, and can 
be written as 

fix^yilx^yj) = ^[A(xi - Xj) - A( Xi - y,) - A( Vi - Xj ) + A( Vl - yj )f (3.2) 

where A(r) is the Green's function given by 

., , l <Pk e' kr , , 

AW = iW^- (33) 

To take confinement into account we replace the infinite range Coulomb potential with a 
screened Yukawa potential. This implies that the Coulomb propagator 1/fc in eq. fl3.3| ) 
is replaced by l/(fe + M 2 ), where M = l/r max is the confinement scale. As a result, the 



four functions A in eq. (3.2) will be replaced by 

with Kq a modified Bessel function. For small separations, where r <C r max , the function 
-Ko( r / r max) behaves like ln(r max /V), and we then recognize the result in eq. fl2,2|). For large 
separations, r S> r max , Ko(r/r max ) falls off exponentially ~ ^/ " r ™ e" r ''' m ", as expected 
from confinement. 

In a similar way, the underlying Coulomb potential in the dipole splitting function 



in eq. (2.1) can be replaced by a screened Yukawa potential, using again the replacement 
1/k — ► l/(fc + l/r max ). The modified splitting probability is then given by 

^^^ d2 JJ_x^ K f\x-z[\_J_y-z Ki f\y 



CLi Z7T ^ ''max |«" ^>| \ ''max / '"max 12/ 

(3.5) 

For small arguments Ki(r/r max ) ~ r mSLX /r, from which we get back the result in eq. ([2.1|), 
while for large arguments K\(r/r max ) ~ \/^r max /r • e~ r ' rmax , and once again we obtain 
an exponentially decaying field. 

4. Initial wave functions 

4.1 Proton wave function 

In ref. [2] we also introduced a simple model for the proton in terms of three dipoles 
with extensions determined by a Gaussian distribution. The resulting model was in good 
agreement with total cross sections for both DIS and pp collisions. It was shown in [15] 
that the three valence quarks in a proton emits gluons in transverse space with the same 
distribution as three dipoles, only with half the intensity. Thus, by modeling the proton 
with a closed chain of three gluons we emulate the fact that a proton at rest may contain 
more charges than its valence quarks. This is analogous to the finding of Gliick, Reya and 
Vogt, who needed a large valence-like gluon component when trying to fit parton densities 
evolved from a very low scale [16]. Thus, although not a fully realistic description of the 
initial non-perturbative proton state, the model appears to give a fair representation of the 
multi-dipole system obtained at the low x-values, which are important for the high energy 
collisions. 

The results turned out to be almost independent of the shape of the three starting 
gluons, except for the size of the triangle. In fact, equilateral triangles that were allowed to 
vary in only size and orientation turned out to model the proton as well as more complicated 
formations. With a Gaussian distribution for the size of equilateral triangles, motivated by 
the exponential dependence on t for the elastic cross section, data on total cross sections 
for DIS and pp collisions are well reproduced, when the width of the Gaussian was tuned 
to 3.5 GeV" 1 « 0.66 fm. 

As discussed above, the differential and elastic cross sections are determined by the 
fluctuations in the scattering amplitude, and in ref. [3] it was pointed out that a Gaussian 



wavefunction as discussed above must overestimate the fluctuations of the incoming state in 
its rest frame. The probability for the three quarks to simultaneously be located in a single 
point ought to be suppressed, and it was emphasized that the exponential ^-dependence 
of the elastic cross section, which motivated the Gaussian shape, is only observed for 
\t\ < 0.15 GeV , corresponding to b ^ 1 fm. A wavefunction of the form 

\^ = Ce- (r - R " )2/w2 (4.1) 

was also tested, and found to give essentially identical total cross sections. The fluctuations 
are here suppressed by a small value of w, and in ref. [3] it was observed that reducing the 
fluctuations to a minimum gave good agreement with the integrated elastic and diffractive 
cross sections in pp collisions. Lacking further constraints we could, however, only present 
an upper limit for a e \, by neglecting the fluctuations in the wavefunction, thus including 
only those in the cascade evolution. 

An essential motivation for the present analysis of quasi-elastic j*p cross sections and 
of the t-dependence of the pp elastic cross section, is to check whether the fluctuations in the 
dipole cascade model are also consistent with these observables, and if more constraints can 
be put on the shape of the initial proton state. In this analysis we will use the two-parameter 



form in eq. (4.1), and see if this can be adjusted to reproduce also the (quasi-)elastic cross 
sections. 

At this point we also note that in many analyses the fluctuations in the cascade evo- 
lution are neglected. This means that e is replaced by e~ in eqs. (2^) and (|2.5|). 



Including also the fluctuations in the cascade implies that the impact parameter profile has 
to be more "gray" and less "black and white". As an example the amplitude {T(b = 0)) 
is a factor 2/3 smaller in our formalism than in the analysis by Kowalski and Teaney [17], 
for a dipole of size 2 GeV -1 and x = 10 -4 — 10 -5 . 

4.2 Photon wavefunction 

4.2.1 Large Q 2 

For large Q 2 the coupling of the 7* to the qq pair can be calculated perturbatively. The 
well known result to leading order is 



^f h - h (Q,r,z) = v ™ ° Qz(l-z)e f K (re f )8 h - h 



^(Q,r,z) =, y/aEM J Tc, " < f (ir" , [:A h _A li _ - , I - :)4-<U) ~7^(r: 7 ) ! ( 1.2) 



IT 

+5 h+ 5f l+ mfK 

with 

E f = ^z(l-z)Q 2 + m 2 . (4.3) 

Here A = and + denote the longitudinal and transverse wavefunctions respectively, / 
denotes the quark flavour, and Kq and K\ are modified Bessel functions, ej is the electric 
charge of the quark in units of the proton charge and rrif the effective mass of the quark. 



4.2.2 Smaller Q 2 

For smaller Q 2 the photon has also a hadronic component. In [18] it was shown that also 
the total 7*p cross section at HERA could be well described over a wide range of ener- 
gies and virtualities, when the hadronic component was simulated by a relatively small 
effective quark mass ~ 60 MeV. For the exclusive reactions studied here we need a more 
careful treatment of the hadronic component, and we expect that these processes can pro- 
vide relevant constraints on the photon wavefunction. The hadronic component should be 
particularly important for the real photons produced in Deeply virtual Compton scattering 
(DVCS). 

For small Q 2 a small effective quark mass allows for rather large dipoles with a corre- 
sponding large cross section. In the present analysis we include an improved description 



of confinement effects in the dipole evolution (see section 3J3), and we will therefore try to 
include confinement effects also in the photon wavefunction. Our photon model is inspired 
by the Vector Meson Dominance modeling introduced by Forshaw et al. in [9] (which, in 
turn, was inspired by [19]). This model contains an enhancement factor for dipoles of a 
typical hadronic size, together with a large quark mass which suppresses dipoles larger 
than the confinement scale. In our model we will use the same enhancement factor, but 
we use a suppression of large dipoles related to the confinement scale r max , instead of the 
large quark mass used in ref. [9]. 

The actual implementation in our MC program relies on shrinking dipoles larger than 
'"max by reducing the size r pert generated according to the perturbative photon wavefunction 
to r soft , defined by 



'"soft \f port J — -"shrink 



N 



\ -""shrink / 



For small dipoles this gives r so f t ss r pert , but for large dipoles it gives a Gaussian suppression. 
The parameter -R sn rink is adjusted to give the same effective cutoff, r max , as the one obtained 
for large dipoles in the cascade evolution. This is obtained for -R sn rink = 4.3 GeV. 

The enhancement factor for dipoles with a typical hadronic size, introduced in ref. [9], 
is given by the form 

f(r ,_ l + B v eM-{r-Rv?/w 2 v ) 

J[) l + B v exp(-R 2 v /w 2 v ) { ' 

This factor multiplies the squared photon wavefunction after the shift in eq. (fOI). The 
enhancement resembles very much the shape we use for the proton wavefunction in eq. (|4.1f ) , 
and we can think of the the whole correction 

Import)! 2 -► |* 7 (r so ft)| 2 /(r so ft) (4.6) 

as the modeling of the virtual photon fluctuating into vector meson states with r-values of 
a hadronic scale. Partly this enhancement can be thought of as due to a longer lifetime of 
these states, and partly a simulation of a gluonic component in the vector meson, in a way 



similar to our model of the initial proton wavefunction in section 4.1 . The photon model 
contains three adjustable parameters, By, Ry, an d wy, which have to be determined from 
experiments. 
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Table 1: The parameters used for the boosted Gaussian and DGKP wavcfunctions in this paper 
in GeV-based units. 

4.3 Meson wavef unctions 

The wavefunction of a vector meson cannot be calculated perturbatively, and has to be 
described by models. In the rest frame it is generally assumed that the lowest Fock state 
with a single qq pair dominates. This component must then be normalized to 1, in con- 
trast to the photon for which the qq state is a perturbative fluctuation. In addition the 
wavefunction at the origin is determined by the decay rate of the vector meson. Thus there 
are two constraints allowing two parameters in an ansatz to be determined. In a boosted 
frame higher Fock states may then be generated by gluon emission. Different models can 
differ in the functional form used for the wavefunction in the rest frame, and in the de- 
scription of the transition from coordinate space to the momentum fractions z and 1 — z 
used in the light-cone wavefunction. We will here concentrate on the DGKP model [20] 
and the boosted Gaussian model [9], which in the analysis by Forshaw et al. give the best 
agreement with the experimental data. 

4.3.1 The DGKP model 

In this model for the meson wavefunction, proposed by Dosch, Gousset, Kulzinger, and 
Pirner [20], it is assumed that the dependence on the transverse and longitudinal coordi- 
nates, r and z, factorizes. The transverse part of the wavefunction is assumed to be a pure 
Gaussian, consistent with soft hadron-hadron collisions. For the longitudinal component 
it assumes the form proposed by Wirbel, Stech and Bauer [21]. The resulting light-cone 
wavefunction has the following form: 



*/&(»-.*) 

*/£('.*) 



,^My5_ h - h z{l - z)y/z(l-z) 



exp 



r 2 uj 2 T 



exp 



M v {z - 0.5) 5 



.yV+ (u^rie i6 \z5 h+ 5 h _ - (1 - z)5 h -5- h+ ) + mA yj z{l - z) x 
M 2 (z -0.5) 2 ' 



exp 



2u;f. 



exp 



r 2 io% 



(4.7) 



Here My is the mass of the vector meson, and the size parameter uj and the normalization 
constant JV are determined from the electronic decay rate and the normalization condition. 
(Our notation differs from the original paper, as we have collected the multiplicative factors 
in the normalization constant r yV.) The shape of the wavefunction of the p with the 
parameters we have used (see table [l]) can be seen in figure |]. 
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Figure 1: The wavefunctions l^^r, z)\ (left) and |^t(^, z )\ (right) of the DGKP model for a 
p meson with our quark mass of 60 MeV. Note the different scales in r, both in GeV -1 . 

4.3.2 The boosted Gaussian model 

The "boosted" models are obtained by assuming a given wavefunction in the meson rest 
frame. This is then boosted into a light-cone wavefunction using the Brodsky-Huang- 
Lepage prescription [22], in which the invariant mass of the quark-antiquark pair is the 
same in the rest frame and the light-cone frames. The result of this procedure is not 
factorizing in r and z. In the simplest version the initial wavefunction in the rest frame 
is a simple Gaussian. In an alternative version by Nemchik et al. [23] a hard Coulomb 
contribution is added, dominating for small r. For the pure Gaussian version suggested 
by Forshaw et al., which we assume in this analysis, the resulting wavefunction has the 
following form 



*tk(r>*) 



fhh 



**&<r>*) 




z)My + m 2 



Xi 



f 



R? 



mJR 2 
z(l-z) 



cxp 






-2z(l -z)—z exp — 



Z )^2 ie * ( z5 h+k- 



R 2 J 

[1 - z)5 h -5- h 



z(l — z)r s 

R o 

m f 



'hi, 



(4.8) 



m)R 2 



exp I —2z(l 



'R? 



cxp 



) + mfSh-\ 
R 2 



-m) 



In this model the transverse size R of the meson, and the normalization jV\ are the two 
parameters to be tuned. (In our notation all multiplicative constants have also here been 
included in a single normalization constant.) The shape of the wavefunction for the p 
meson using the parameters in table is shown in figure g. 

5. Tuning of parameters and the differential pp cross section 

5.1 The total and elastic pp cross section 

We start by tuning the model to pp scattering data. Here the model contains 4 parameters, 
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Figure 2: The wavefunctions |\Pl(»", z)\ (left) and \^T(r, z)\ (right) of the Boosted Gaussian 
model for a p meson with our quark mass of 60 MeV. Note the different scales in r, both in GeV -1 . 



Aqcd and r max describing the dipole evolution, and R p and w determining the proton wave 
function in eq. ( f4,ip (with C fixed by normalization). In ref. [3] we found that the values 
for Aqcd and r max are correlated, such that a larger r max can be compensated by a smaller 
Aqcd- It was also noted that the integrated elastic cross section favors a narrow proton 
wave function, corresponding to a small value for the parameter w. A large upvalue, or 
a single Gaussian oc exp(— r 2 /R 2 ), gives too large fluctuations and correspondingly a too 
small elastic cross section. To constrain the fit we here add the differential cross section 
da/dt to the integrated cross sections <7 t ot, o"diff> and a c \ studied in [3]. We will then in 
next section check if the result also can reproduce the quasi-elastic cross sections in j*p 
collisions. 

With the proton wavefunction given by eq. (|4.1|) the total and elastic cross sections 
are given by 



dtot = 2 / d 2 bd 2 r pl d 2 r p2 \^ p (r pl )\ 2 \^ p (r p2 



)| 2 (1 



112, 



(5.1) 



Ccl 



d 2 b 



d 2 r n id 2 r„ 2 \y 



pia r p2 I 



<"pU 



l*p(r p2 )| 2 (l 



e- F )i 2 



(5.2) 



Here b is the impact parameter, v p i (i = 1,2) parameterizes the size and orientation of 
the triangles describing the colliding protons. The Monte Carlo is used to simulate the 
dipole evolution in the rest frame of the collision, and to calculate 1 — e~ F . The average 
(1 — e _i? )i2 is over simulations for fixed impact parameter and starting dipole states ri 
and r 2 . Note that in the elastic cross section the average over evolutions and the integrals 
over the wave functions is taken on amplitude level before taking the absolute square. 

When tuning the parameters we find that all observables are almost independent of 
w below 0.5 GeV" 1 . We therefore decided to neglect the fluctuations in the proton size 
completely and set the width to zero, turning the proton wavefunction into a delta function 

at lXr\- 
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Figure 3: The total and elastic cross section for pp collision. The upper cross sections are total 
cross sections, while the lower cross sections are the elastic ones. Tevatron data are from [24-27], 
SPS data are from [28] and cosmic ray data are from [29]. The lines are our model with tuned 
parameters. 

If the total and integrated elastic cross sections are tuned at one energy, we find that the 
energy dependence of these cross sections depends very weakly, if at all, on the parameters 
of the model. Thus this energy dependence cannot be tuned, and the fact that it is close 
to the experimental results is therefore a direct consequence of the model. Our results for 
the total and elastic pp cross sections can be seen in figure y. 

Extrapolating to higher energies we find the total cross section at the LHC nominal 
energy, 14 TeV, to be about 125 mb (117 mb at 10 TeV). We note that this is a rather 
high value compared to other predictions. Thus the Donnachie-Landshoff parameterization 
gives 101.5 mb at 14 TeV [30], while an analysis by Khoze, Martin, and Ryskin gives about 
90 mb [31]. The predicted elastic cross section is about 31 mb for the LHC at 14 TeV 
(28 mb at 10 TeV). 

5.2 The differential elastic pp cross section 

The differential elastic pp cross section is given by 



dt 



1 

4-7T 



d 2 he 



***(£ 



r pl d 2 rp 2 \^p(r pl ) 



l* P (r p2 )| 2 (l 



e~ F h 2 



with t 



(5.3) 

We here neglect the real part of the amplitude, and therefore da /dt will have zeroes from the 
Fourier transform of the amplitude in eq. ( |5.3j ) . Even though the true complex amplitude 
will not be identically zero, the real part is still assumed to be small, producing a dip at 
some value t = to, related to the inverse square of the size of the proton at the relevant 
energy. This dip is visible in some of the experimental data shown in figure [4], where 
we have also included the results from the simulations. The parameters which are most 
sensitive to these distributions are R p , which determines the size of the proton at rest, 



13 




-t (GeV^ 



Figure 4: The elastic cross section as function of t in mb/GeV 2 . The numbers in parenthesis 
indicate how the data has been scaled. The lines are our model with tuned parameters. Predictions 
for 14 TeV is also included. Data arc from [24,32], [33] and [34]. 



and r max which regulates the maximal size of new dipoles, and therefore the increase with 
energy of the proton radius and the variation of the dip position. However, the slope of 
the distribution is basically independent of our parameters, as is the cross section at large 
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Aqcd 


0.2 GeV 


'"max 


2.9 GeV" 1 


R p 


3.0 GeV" 1 


w 


GeV" 1 



Table 2: The tuned values the parameters for the proton wavefunction and the perturbative 
evolution used for our model throughout this paper. 

t-values. Nevertheless, we are able to get a very good description of the data at all i-values 
even though the cross sections vary over many orders of magnitude. In figure || we also 
show our result for the LHC, which predicts the location of the dip in the ^-dependence at 
0.43 GeV 2 at y/s = 14 TeV (0.47 GeV 2 at 10 TeV). The values of the tuned parameters 
can be found in table 0. 



5.3 The total 7*p cross section and tuning the photon wave function 



We will here use ^<y(Q, r, z) to denote the photon wavefunction in eq. ( ^4.6[) , where for small 
Q 2 the perturbative wavefunction is modified to account for the hadronic component of 
the photon. The total j*p cross section can be written 



vtot(Tf*p) 



Xf J 



d 2 hd 2 r p d 2 rdz \^ p (r p ) 



%j{Q,r,z) (1 



6 /dpi 



(5.4) 



where A is the polarization of the photon and / is the flavour of the quark-antiquark pair 
created by the photon. (1 — e~ F )d p is now an average of the evolution of the dipole from 
the photon side and of the dipoles from the proton side. It depends not only on the total 
energy, the size of the proton and photon dipoles and b, but has also a weak dependence 
on z. 

The three parameters By , Ry and wy in the enhancement factor in eq. ( |4.5| ) were fitted 
to the total 7*p cross section as measured at HERA. Here the value of Ry determines the 
range in Q 2 where the enhancement is significant, while wy determines how fast it falls off 
for large Q 2 . The parameter By is just an overall strength of the hadronic component. 

A good fit to data was obtained with a wave function for the hadronic component 
similar to the proton wave function, having a size Ry ~ 3 GeV -1 and a small width. 
The total 7*p cross section with and without the effects of confinement and vector meson 
contributions are shown in figure || The tuned values are given in table ||. These parameters 
are quite different from the ones used by Forshaw et al., who had an Ry of 6.84 GeV -1 
and a fairly large width [35], which thus gives a stronger enhancement for large dipoles. 
For large dipole sizes the elementary dipole-proton cross section assumed in ref. [9] are 
also significantly different from the corresponding ones in our model. The reason why we 
anyhow can get similar results is that in [9] the very large dipoles are suppressed by a 
large quark mass, and the enhancement therefore not very effective. A reason for a smaller 
width in our wavefunction is also that fluctuations in the dipole cascade are included in 
our formalism, which is compensated by less fluctuations in the wavefunction. 
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Figure 5: The total cross section of photon-proton collision as function of photon virtuality and 
center-of-mass energy. The dashed line is calculated with a purely perturbative photon wavefunc- 
tion, while the full line is with a photon wavefunction with both confinement and VMD corrections. 
Experimental data are from [36] 



R v 


3.0 GeV" 1 


wy 


0.2 GeV" 1 


Bv 


9.0 



Table 3: The tuned values of the parameters of the vector meson resonance function /(r) used for 
our model throughout this paper. 

We also note that the cross sections presented in fig. [| are somewhat smaller than the 
corresponding results in refs. [3] and [18]. This is a consequence of the more consistent 
treatment of confinement in the present analysis, which gives a stronger suppression for 
larger dipoles. We believe, however, that a much stronger test of the hadronic component 
will rely on the results for the quasi-elastic reactions discussed in next section. 



6. Results for quasi-elastic 7*p collisions 

In this section we will study predictions for quasi-elastic "y*p collisions, using the photon 



wavefunction parameters determined in section 5.3 



6.1 Deeply Virtual Compton Scattering 

In Deeply Virtual Compton Scattering, DVCS, the incoming particle is a virtual photon 

7./ 



while the outgoing particle is a real photon with wavefunction $;* JQ,r,z). According to 



eqs. (2J3, 2J}) the cross section is given by 



odvcs 






d 2 r p d 2 rdzY,\^p(.r p )\ 2 ^ f ' rX (Q,r,z)^}\0,v,z)(l-e- F } dp 
f 



(6.1) 
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Figure 6: The cross section of j*p — ► 7p for W — 82 GeV as function of Q 2 (a) and as function 
of W for three different Q 2 (b). Experimental data are from [37,38]. 
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Figure 7: The cross section of 7*p — > 7J> as function of t. The brackets refers to (W, Q 2 )- Note 
that the three series have been scaled by 10,1 and 0.1 for better readability. Experimental data are 
from [38]. 



The differential cross section da/dt is obtained from the Fourier transform as shown in 
eq. ( |2.10|) . The results obtained using the parameter values in table || are presented in 
figures y and 0. We see that the results from the model indeed agrees with the data, both 
in normalization and in the dependence on Q 2 , W and t. As this quasi-elastic reaction is 
very sensitive to the fluctuations and the transverse size of the hadronic component of the 
real photon, this is a clear support for the proton-like wave function. 



6.2 Exclusive Production of Light vector Mesons 

The cross section for exclusive vector meson production, rfp — > Vp, can be calculated in 
exactly the same way as for DVCS, only replacing the real photon wavefunction by a meson 
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Figure 8: The cross section for j*p — ► pp. (a) As function of the photon virtuality. The HI data 
has been moved down a factor 10 for better readability. The full and dashed line are with the two 
different meson wavefunction described in |4.3| . (b) As function of the center-of-mass energy W. 
Experimental data are from [39-41]. 

wavefunction: 



Ccl 



/ >b]T f dhd\dzY,\%(r p )\ 2 V?\Q 2 ,r,z)yy\r,z)(l- e- F ) dp 
J x J f 



. (6.2) 



As before we have ignored the real part of the amplitude. Contrary to the case of pp 
scattering it has been shown in [9] that in exclusive production of light vector mesons the 
real part can be quite large, for large Q 2 or large W as large as one half of the imaginary 
part. This would mean that we underestimate the cross section by up to 25% in these 
regions. However, compared to the uncertainties in the vector meson wavefunctions, this 
error is small. 

From our tuning of the hadronic part of the photon wavefunction, it could seem nat- 
ural to assume some universal hadron size and maybe try to model the vector meson 
wavefunctions as a simple gluon-gluon dipole with a size of 3 GeV -1 and a small width. 
However, this would not naturally give us a z-dependence and we would not include the 
possibility that the vector meson may fluctuate into a photon, which could correspond to 
an enhancement at small r-values. Therefore we will simply use the boosted Gaussian and 
DGKP wavefunctions introduced in section 4.3 to estimate the "f*p — > Vp cross section. 



Throughout we will use the parameters listed in table [l]. 

As before, the t-dependence of the cross section can be calculated through a Fourier 
transform of the amplitude. We are also able to calculate the ratio between the longitudinal 
to the transverse cross sections and compare with experimental data. 

Starting with p meson production, the results are shown in figures |8| and |9[ We see that 
the model calculations reproduce experimental data rather well, including the dependence 
on virtuality Q 2 , energy W and momentum transfer t. The Boosted Gaussian wavefunction 
is in general providing the closer fit, while DGKP is having problems mainly in the ratio 
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Figure 9: (a) The differential cross section for -f*p — > pp as function of transferred momentum 
squared \t\. Three different Q 2 has been used, all with a center-of-mass energy of 90 GeV. (b) The 
ratio of longitudinal and transverse cross section for j*p — > pp as function of the photon virtuality. 
Experimental data are from [39,40]. 
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Figure 10: The cross section for j*p — > <j>p (a) As function of the photon virtuality at center- 
of-mass energy 75 GeV. The full and dashed line are with the two different meson wavefunction 
described in section |4.3[ (b) As function of the center-of-mass energy for four different photon 
virtualities Q 2 . Experimental data are from [42,43]. 



between longitudinal and transverse cross sections. It should be noted, however, that this 
ratio could be changed if we decided to use different parameters for the resonance function 
in eq. ( |4.5| ) for the different polarizations. 

Also in (p production our model agrees well with experimental data, as can be seen in 



figures |10| and 11, The qualitative properties are similar to those of p production 

We note that the more stringent test of the hadronic component of the photon is 
obtained from DVCS. The ratio between vector meson production and DVCS is then more 
a test of the vector meson wavefunctions. It is therefore not surprising that we here get 
results similar to those in ref. [9]. The ^-dependence presented in fig. || is, however, a result 
which in our model is sensitive to both the photon and vector meson wavefunctions, while 
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Figure 11: The ratio of longitudinal and transverse cross section for j*p — ► <frp as function of the 
photon virtuality at a center-of-mass energy of 75 GeV. Experimental data are from [42,43]. 

in ref. [9] it was fixed by experimental data. From fig. ^ we see that for lower Q 2 the 
slope in the model is somewhat too steep, thus indicating a too wide wavefunction for the 
p meson. We see from figs, [I] and [2] that the p wavefunctions for transverse polarization 
are extending well beyond 5 GeV - . A faster falloff for large r-values would here give a 
better agreement with the observed t-dependence. 

6.3 Exclusive tp Production 

In the case of i/j production we necessarily encounter more uncertainties. The result is 
sensitive to the mass of the charm quark, and here we use the value 1.4 GeV, which in 
the analysis in ref. [18] gave the correct charm contribution, i 7 ^, to the proton structure 
function. We have not included a ip component in the photon wave function, although this 
would in principle be possible. For the if) meson wavefunction we use the parameters in 
table |. 

Our results are presented in figures |l^ (dependence on Q 2 and W) and |l3| (dependence 
on t). We note that the Boosted Gaussian wavefunction gives a too low cross section over 
the entire kinematic region, while DGKP model agrees very well both in its normalization 
and in its dependence on Q 2 and W. Both models show, however, a somewhat steeper 
t-dependence than the experimental data. This can be compared to the corresponding 
result for p meson production, and we conclude that also for the ip meson the parameters 
in table [l] gives somewhat too wide wavefunctions. 

7. Conclusions and Outlook 

In this paper we have spent some effort on the modeling of non-perturbative aspects of 
the proton, photon and vector meson wavefunctions. None of our models are in any way 
unique or on solid theoretical grounds. However, they do allow us to compare our dipole 
evolution model directly to experimental data. Fixing the wavefunction parameters at one 
energy we find that the energy dependence of both total and (quasi-) elastic cross sections 
are well described by the cascade evolution, and rather independent of our modeling of the 
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Figure 12: The cross section for 7*p — v tpp (a) As function of the photon virtuality at center- 
of-mass energy 90 GeV. The full and dashed line are with the two different meson wavefunction 
described in |4.3| . (b) As function of center-of-mass energy for Q 2 = 22 A GeV 2 (lower data) and 
3.2 GeV 2 (upper data). Experimental data are from [44]. 



> 

CD 

o 

!3 



T 

I 

o 

"D 



0.2 





\ 1 1 1 1 

H1 


o 


100 


< \4> BG - 


- 




: \ \ DGKP - 


■ 




: \"H 


- 




\. \ 






K \\ 




10 


\K \\ + 


~ 




: \/K ^\\ 


" 




^\ ""- ^\. 


^"^ 


-1 


I \l H. T | 


^ 



0.4 
-t (GeV 



0.6 

2> 



0.8 



Figure 13: The cross section for j*p — > ipp as function of transferred momentum squared t for 
Q 2 = 22 A GeV 2 (lower data) and 3.2 GeV 2 (upper data). Experimental data are from [44]. 



wavefunctions. Also the slope in da/dt for elastic pp scattering and DVCS is in agreement 
with experimental data independently of the tuning. This indicates a very high predictive 
power of our evolution model both when it comes to the average multiplicity and sizes of 
the dipoles (mainly influencing the total cross sections), the rate of increasing transverse 
size due to the dipole cascade (determining the energy variation of the dip position), and 
the fluctuations (mainly influencing the magnitude of elastic cross sections and their t- 
dependence). 

Nevertheless, our modeling of the non-perturbative wavefunctions does give us valuable 
insights. Including the fluctuations in the cascade, the fluctuations in the proton wave- 
function have to be rather small, in order to give the observed elastic cross section. The 
photon wavefunction clearly needs a hadronic component with a wavefunction with similar 
size as the proton and with similarly small fluctuations. The fact that the size comes out 



21 



to be the same as the size of our proton may be a coincidence, but it could also indicate 
that there is a universal size of hadrons, at least when consisting of light quark flavours. 

For the vector meson wavefunctions we have tested two different forms, which in the 
analysis by Forshaw et al. gave the best agreement with data for diffractive vector meson 
production. For light mesons the best result was obtained by the boosted Gaussian wave- 
functions, while for ip production the DGKP wavefunction was favored. In both cases the 
t-dependence was somewhat too steep, indicating that these wavefunctions extend out to 
too large r-values, where in particular the wavefunctions for transverse polarization are 
much wider than our wavefunctions for the proton and the photon. 

The robustness of our model for dipole evolution, both for inclusive and exclusive 
cross sections increases our confidence that it can be used to also model fully exclusive 
final states. In future publications we will therefore concentrate on turning our Monte 
Carlo simulation program into a full-fledged event generator, which would then be able to 
model multi-particle production at e.g. the LHC, with special emphasis on the underlying 
event and multiple scatterings. 
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